Extreme nonlinear electrodynamics in metamaterials with very small linear dielectric 
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We consider a sub-wavelength periodic layered medium whose slabs are filled by arbitrary linear 
metamaterials and standard nonlinear Kerr media and we show that the homogenized medium 
behaves as a Kerr medium whose parameters can assume values not available in standard materials. 
Exploiting such a parameter availability, we focus on the situation where the linear relative dielectric 
permittivity is very small thus allowing the observation of the extreme nonlinear regime where the 
nonlinear polarization is comparable with or even greater than the linear part of the overall dielectric 
response. The behavior of the electromagnetic field in the extreme nonlinear regime is very peculiar 
and characterized by novel features as, for example, the transverse power flow reversing. In order to 
probe the novel regime, we consider a class of fields (transverse magnetic nonlinear guided waves) 
admitting full analytical description and we show that these waves are allowed to propagate even 
in media with e < and /i > since the nonlinear polarization produces a positive overall effective 
permittivity. The considered nonlinear waves exhibit, in addition to the mentioned features, a 
number of interesting properties like hyper-focusing induced by the phase difference between the 
field components. 

PACS numbers: 78.67.Pt, 42.65.Tg 



I. INTRODUCTION 



Electromagnetic propagation through metamaterials 
has stimulated, in the last decade, an intense research ac- 
tivity with two main purposes: the identification of artifi- 
cial structures exhibiting anomalous values of the permit- 
tivity e and the permeability [i [HQ and the investigation 
of the electromagnetic phenomenology resulting from the 
unusual electromagnetic properties. Leading examples of 
the novel class of phenomena supported by metamate- 
rials are superlensing [3, optical cloaking [j| 0| and 
photonic circuits @- Recently, a good deal of attention 
has been devoted to the investigation of metamaterials 
characterized by a very small dielectric permittivity (e- 
near-zero materials) since they host an electromagnetic 
regime where the magnetic field displays static features 
Q and the tunneling of electromagnetic energy through 
sub-wavelength channels has been proposed [9( and ex- 
perimentally observed [10|, . Metamaterials exhibiting 
remarkable nonlinear properties have been investigated 
as well [l2| and various soliton manifestations [13-[l6( 
have been considered. The investigation of metamaterial 
nonlinear properties is particularly important in that it 
can lead to overcoming one of the fundamental limit of 
nonlinear optics, the fact that most of the optical materi- 
als have a relatively weak nonlinear response. The main 
idea is that the local electromagnetic fields of the inclu- 
sions in the metamaterial can be much larger than the 
average value of the field thus producing an enhancement 
of the nonlinear response [17Hl9l | . The problem of achiev- 
ing a substantial enhancement of the nonlinear response 
has also been considered within the more general subject 



of composite structures homogenization [2CH24{ , and the 
strategy is always that of conceiving a microscopic inho- 
mogeneous structure concentrating the field within the 
nonlinear constituents. 

The full exploitation of the nonlinear response is pos- 
sible only if the nonlinear polarization is not a small per- 
turbation to the linear part of the electric displacement 
field and generally this is achieved through nonlincarity 
enhancement or by means of resonant processes or pho- 
torefractive processes where the large nonlinearities come 
at the cost of a large time response. However, as shown 
in this paper, the interplay between the linear part of 
the electric displacement field and the nonlinear polar- 
ization can be made efficient even by following the oppo- 
site route, i.e. by reducing the linear polarization. We 
therefore devise a nonlinear medium with a very small 
linear dielectric constant since it is a natural setting for 
the observation of the electromagnetic regime where the 
nonlinear response does not play the role of a mere per- 
turbation. 

In this paper we consider a periodic layered compos- 
ite whose slabs are filled either with linear media with 
arbitrary permittivity and permeability or by standard 
isotropic focusing or defocusing Kerr media. Exploiting 
a suitable extension of the well known technique generally 
used for describing the homogenization of linear layered 
composites, we show that the homogenized medium is 
characterized by effective constitutive relations formally 
coinciding with those of a standard Kerr medium. We 
note that the parameters characterizing such an effec- 
tive response can be simply tailored through a suitable 
choice of the composite underlying constituents and that, 
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due the large freedom in choosing both constituent me- 
dia and their volume filling fraction, the design of lin- 
ear and nonlinear properties can be independently per- 
formed. Therefore our composite medium allows a full 
and efficient engineering of the Kerr nonlinear response. 
More interestingly, we prove that the effective response 
parameters span very wide ranges encompassing values 
not available in standard media. The standard isotropic 
Kerr response (in the the frequency domain) generally de- 
pends on two parameters usually denoted with A and B 
[26j in term of which we define \ — A and 7 = B /2A and 
the available values of 7 belong to the range < 7 < 3 
(depending on the actual mechanism supporting the Kerr 
response). In this paper we prove that the parameter 7 
appearing in the effective Kerr nonlinear response can, in 
principle, assume any value (encompassing negative and 
very large values) and we discuss the impact of the exotic 
values of 7 on the electromagnetic phenomenology. 

Confining our attention to transverse magnetic (TM) 
fields, we show that it is possible to design a nonlinear 
Kerr medium whose linear dielectric permittivity is much 
smaller than one and therefore able to host the electro- 
magnetic regime where nonlinearity can not be regarded 
as a perturbation (extreme nonlinear regime). As a first 
general electromagnetic effect characteristic of the ex- 
treme nonlinear regime, we discuss the transverse power 
flow reversing, i.e. the fact that, for an electromagnetic 
beam, the Poynting vector on the beam propagation axis 
can be antiparallel to Poynting vector on the beam lateral 
sides [25|]. Within the extreme nonlinear regime, the full 
exploitation of the nonlinear response can be achieved 
and, in order to discuss the consequent novel electro- 
magnetic phenomenology, we consider a class of nonlin- 
ear guided waves admitting full analytical treatment. A 
number of the obtained nonlinear guided waves exhibit, 
as expected, transverse power flow reversing. We prove 
that, in the situation where e < 0, ^ > and \ > (or 
in the equivalent case e > 0, [i < and x < 0) these 
waves are allowed to propagate through the medium if 
— 1 < 7 < 0, in striking contrast to the fact that, in a lin- 
ear medium with the same permittivity and permeability, 
any wave would be evanescent. This is possible since non- 
linearity overcomes the linear contribution to the electric 
displacement field producing an effective nonlinear di- 
electric permittivity able to support propagating waves. 
In addition we identify a mechanism which, combining 
the extreme nonlinear regime with the properties of the 
effective medium response with 7 ~ —1, produces an 
hyper-focusing of the considered nonlinear guided waves. 
In the situation e > 0, /i > and x < (or in the 
equivalent case e < 0, fi < and x > 0) the nonlinear 
guided waves exist in the wider range — 00 < 7 < — 1 and 
< 7 < +00 thus allowing to probe the Kerr nonlinear 
response for large values of |-y|. Since the phase difference 
between the components of the considered waves is 7r/2, 
we predict that, for 7 — > +00, an extreme compensation 
occurs (within the effective nonlinear polarization) and 
supports a field whose transverse E x and longitudinal 
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FIG. 1: Geometry of the nonlinear layered composite and TM 
electromagnetic field configuration. 



E z components are such that \E X \ ~ \E Z \, i.e they share 
the same profile. 

The paper is organized as follows. In Sec. II we inves- 
tigate the homogenization of a one-dimensional periodic 
layered medium comprising both linear (metamaterials) 
and nonlinear (Kerr media) slabs and we derive the con- 
stitutive relations characterizing the effective medium. In 
Sec. Ill we focus on the transverse magnetic field config- 
uration and we discuss the effective nonlinear response 
engineering and the allowed extreme nonlinear regime, 
together with the transverse power flow reversing effect. 
In Sec. IV we investigate a class of nonlinear guided waves 
belonging to the extreme nonlinear regime and we discuss 
a number of their peculiar properties. The existence of 
the considered nonlinear guided waves is investigated in 
Appendix A. 



II. HOMOGENIZATION OF A 1-D NONLINEAR 
LAYERED MEDIUM 

Consider a monochromatic electromagnetic field 
(whose time dependence is exp(— iuit)) propagating 
through a metamaterial layered medium consisting of pe- 
riodically repeating, along the y— axis, N layers of dif- 
ferent media of thicknesses dj (j = 1...N), so that the 
structure spatial period is d — YliLj 4? ( se e Fig.l where 
the case N — 4 is depicted). Each one of these iV me- 
dia can be either a linear metameterial with arbitrary 
dielectric and magnetic properties or a cubic standard 
nonlinear medium so that, the electromagnetic response 
of the j— th medium is modelled by the general constitu- 
tive relations 

D, = eocjEj + eottPj ■ E*)E, + 7j (E J ■ E^E*], 
Bj = fiofijUj, (1) 

where eo &ud ^0 are the vacuum permittivity and per- 
meability constants, whereas Ej, Dj, H,-, Bj are the 
complex amplitude of the local electromagnetic field vec- 
tors, Cj and [Lj are the relative permittivity and perme- 
ability of the j— th layer, Xj an d 7j are the standard 
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parameters characterizing the isotropic cubic nonlinear 
response of the j—th layer [26]. Evidently, \j = if 
the j—th layer is filled by a linear medium and fij = 1 
for a nonlinear dielectric. The linear dielectric constants 
€j are here regarded as arbitrary complex numbers since 
the layers can be filled with lossy and active media. If 
the field vacuum wavelength is much greater than the 
spatial period (A = 2irc/u> 3> d) the considered period- 
ically nonlinear stratified medium can be homogenized, 
i.e. its electromagnetic response can be shown to coin- 
cide with that of a suitable homogeneous medium. In 
order to obtain the overall effective response, note that 
any component of the local electromagnetic field can be 
assumed, within each layer, independent on y since the 
layers are extremely small (d -C A). Although this is a 
very reasonable physical assumption, it can be rigorously 
proven exploiting the well known powerful two-scale ex- 
pansion method [13, [Hj]. Each of the physical observ- 
able (macroscopic) electromagnetic field vectors, say V 
(V = E, D,B,H), is obtained by averaging (along the 
y— axis) the layer local fields over the period d, so that 



v v, ^/,v. 



(2) 



where fj — dj/d is the volume filling fraction of the j—th 

medium Q^jLi fj = -0 ana - the averaging has been per- 
formed by exploiting the uniformity along y of the local 
fields Vj . Within each unit cell (of thickness d) , at each 
plane interface between the j—th and (j + 1)— th layer, 
the local fields have to satisfy the electromagnetic bound- 
ary conditions (continuity of the tangential component of 
electric and magnetic fields and continuity of the normal 
components of the displacement and magnetic induction 
fields) so that the local fields are joined by the relations 



Ej x — 


E(j+i) x , 




= H (] + l)x, 




E jz = 


E(j+i)z, 


H jz 


= H {j+l)z, 




Dry = 


D {j+l)y, 


B jy 


= B U+i)y 


(3) 



where j = l...(N — 1). Averaging the local fields Bj X , 
Bj Z and Hj y and exploiting the second of Eqs.(TT]) to- 
gether with magnetic boundary conditions of Eqs.Q, it 
is straightforward to prove that the macroscopic fields 



B 



(Bj) and H 



(H 



are related by 

B = flQUeffH 



(4) 



where the relative magnetic permeability /U e // is the di- 
agonal tensor fj, e ff — diag[(nj), (fij 1 )^ 1 , (fij)], which is a 
very well-established result concerning the homogeniza- 
tion of layered media [Hj]. As far as the dielectric re- 
sponse is concerned, exploiting the first of Eqs.(TT]) to- 
gether with the fact that Ej X — (Ej X ) — E x , Ej Z = 
(Ej Z ) — E z and Dj y = (Dj y ) = D y (obtained by com- 
bining the electric boundary conditions of Eqs.© and 
Eq.(|2J|) we obtain 



D 



)E 



+ e [(x,)(E ± ■ E1)E ± + (x,-7i)(E± ' E±)E±] 

+ ^[{Xj\E jy \ 2 )^ + (XjljE%)^l], 

D y = e ejEj y 

+ e oXj [(E ± ■ El)E 3y + 7i (E ± • V±)E* y ] 



+ eoXj [(l+Jjm^E 



JUi 



(5) 



where D_l = D x e x + D z e z and Ej_ = E x e x + E z e z are 
the transverse parts of the macroscopic displacement and 
electric fields, respectively and D y is the y— component 
of the macroscopic displacement field. Note that in the 
second of Eqs.([5]) no averaging has been performed. As 
a consequence of the layers nonlinear dielectric behavior, 
the first of Eqs.© contains terms where the squares of 
Ej y is suitably averaged. In order to derive the effective 
medium response, the second of Eqs.© has to be solved 
to express Ej y (j = 1...N) as a function of the macro- 
scopic fields Ej_ and D y . This can be pertubatively done 
by noting that, since we are considering the standard 
Kerr effect for each nonlinear layer [26j ]. the terms con- 
taining Xj are much smaller than the term eoCjEj y so 
that, up to the first order in the field cubic terms, from 
the second of Eqs.([5]) we obtain 



e E jy = -D y -^ [(El • El) A, + 7i (E ± • E ± )D* y ] 



Xj- 



■lj 



2 4 

e o e j 



DyYDy. 



(6) 



Averaging Eq.([B]) we obtain a relation joining the macro- 
scopic field Ey and the fields Ej_ and D y so that this 
relation, if its nonlinear contributions are much smaller 
than the leading linear term, can be inverted to pertur- 
batively yield D y as a function of E^ and E y . Therefore, 
up to the first order in the field cubic terms, we obtain 



D y = eoiejY'E, 



(Xjtj 2 ) 



E, -E* 



(X J (1 + 7 J )S T4 ) 

2(^ 1 ) 4 



E„ 



-Ei • E j 



(Xj(l 



E 



E„ 



E* 



(7) 



Substituting the fields Ej y of Eqs.© into the first of 
Eqs.© and using Eq.Q (neglecting everywhere the 
terms containing powers of Xj higher than one) we obtain 



D 



eo(e,-)E_L 



£0 



+ eo 



(Xj)Vx ■ E^ 



(Xjlj)Vi ■ Ej 



E 



EI. (8) 



Equations ([7]) and © solely contain the macroscopic 
fields so that they are the constitutive dielectric rela- 
tions characterizing the effective medium obtained by the 
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homogenization of the considered ID layered medium. 
The obtained dielectric response shows that the effective 
medium behaves like an anisotropic Kerr medium whose 
nonlinear properties can be tailored by suitably choosing 
the underlying layered composite. Note that, in the spe- 
cific case of linear layers (i.e. Xj = for all j), Eqs. ([7]) 
and © yield 



D = e e e// E 



(9) 



where the relative dielectric permittivity e e f / is the diag- 
onal tensor e e ff = diag[(ej) , (ej 1 )^ 1 7 (ej)], reproducing 
a well known result concerning homogenization of linear 
layered media [29|, |3(J ■ 



III. EXTREME ELECTRODYNAMICS OF TM 
FIELDS 

Electromagnetic propagation through the homoge- 
nized medium of Sec. II is described by the macroscopic 
Maxwell equations 



VxE = «wB, 
V x H = -iuT> 



(10) 



and the constitutive relations of Eqs.Q, ([7]) and 
Hereafter, we focus on Transverse Magnetic (TM) elec- 
tromagnetic fields (see Fig.l) of the form 



E = E x (x,z)e x + E z (x,z)e z , 
H = H y (x,z)e y 

for which the constitutive relations become 

D = e eE + e x[(E-E*)E + 7(E-E)E*] J 
B = ^oM H , 

where 



(11) 



(12) 



N 



N 



EM" 1 



3=1 

JV 



(13) 



Therefore the effective dielectric response experienced by 
a TM field coincides with the standard isotropic Kerr 
response. In analogy to what happens in linear lay- 
ered media we note, from Eqs. (ll3[) . that the permittiv- 
ity e and the permeability /j, are the weighted and har- 
monic weighted means, respectively, so that the former 
is bounded by the minimum and maximum of its mi- 
croscopic values whereas the latter is unbounded since 
some fij can be negative. It is worth stressing that the 
averaging of the microscopic linear parameters addition- 
ally allows an efficient loss management since, combining 
lossy and gain media, one can design an effective medium 



whose effective parameters e and \x have negligible imag- 
inary parts HH. On the other hand, x is the weighted 
mean of the constituents nonlinear Kerr coefficients and 
this is in agreement with the results of Ref.fHj]. More 
interestingly, 7 is generally not the weighted mean of its 
microscopic values since some Xj can be negative (corre- 
sponding to defocusing nonlinear layers), so that 7 can 
assume any positive and negative value. This result is 
particularly interesting since, for standard cubic isotropic 
materials, there is, in general, a small number of avail- 
able 7 (depending on the physical mechanism supporting 
the nonlinear response [26|) whereas, for composite ma- 
terials, it has been shown in Ref. [22| that 7 can span the 
whole range < 7 < 3. 

Even though Eqs. (fT2"]) formally coincides with the stan- 
dard isotropic Kerr response, the considered layered 
material can support a nonlinear electromagnetic phe- 
nomenology very different from that observed in a stan- 
dard Kerr medium, as discussed in the following three 
subsections. 



A. Extreme nonlinear regime 

Consider a layered medium for which e is such that 
< Re(e) -C 1, |im(e)| <C Re(e) and the effective non- 
linear coefficient x is °f the same order of magnitude of 
underlying coefficients Xj (i- e - Ixl ~ \Xj\)- An effec- 
tive dielectric permittivity with a very small real part 
and a negligible imaginary part can be achieved by com- 
bining positive and negative standard dielectric layers 
(|i?e(ej)| > 1) together with gain media (generally un- 
avoidable since absorption due to negative dielectrics can 
not be in principle neglected [31(). For the sake of sim- 
plicity let us consider the case 7 = 0. If the electro- 
magnetic field propagating through the medium is such 
that 



\E X 



\E Z 



(14) 



we have |xj(E • E*)| ~ |e| -C \ej\ so that, in the first 
of Eqs.(fT]), the nonlinear part is much smaller than the 
linear contribution eoCjEj. Therefore, the nonlinear lay- 
ers are in the presence of a field whose intensity is suffi- 
ciently small for their response to be purely cubic and, as 
a consequence, the microscopic responses of Eqs.(fT|) hold 
and the macroscopic constitutive relation in the first of 
Eqs. (fT2")) . holds as well. On the other hand, combining 
Eq. (fT4j) and the first of Eq. (fT2"|) we conclude that, as op- 
posed to what happens in standard nonlinear materials, 
the considered layered medium is able to support the ex- 
treme nonlinear regime where the electromagnetic field 
is such that the linear contribution eE and the nonlin- 
ear term x(E • E*)E in the overall dielectric response 
have the same order of magnitude. It is evident that 
the same argument, with some slight changes, allows to 
prove that the extreme nonlinear regime is observable 
for any value of 7. Even though the advantages of such 
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an extreme nonlinear regime are self-evident, it is worth 
to compare it to the standard paraxial nonlinear optical 
situation where the overall refractive index is uq + Sn 
where hq is the background linear refractive index and 
Sn = 712I << no is the Kerr nonlinear term (where I is 
the optical intensity). In this situation the field satisfies 
the equation 



V 2 E + kl{nl + 2n <5n)E = 



(15) 



(fco = oj/c) and, within the paraxial regime, the field 
is of the form E = e lkz A. At this point one chooses 
k = kgno since, substituting the paraxial field into 
Eq.(TT5j) the term A^ngE (responsible for the fast spatial 
variation of the field) is removed and the paraxial equa- 
tion + j^\/ 2 ± A — — ^SnA is readily obtained by ne- 
glecting the term containing d 2 A. Therefore, in standard 
paraxial nonlinear optics, Sn can play an important role 
in Eq. (|15|) since the large linear contribution proportional 
to n 2 , is suppressed by the presence of the carrier plane 
wave. In other words, the nonlinear Kerr contribution 
eoX|E| 2 E is not requested to compete against the whole 
linear part eoeE of the dielectric response. As opposed, 
in the case of the extreme nonlinear regime discussed in 
this section, this competition can happen and therefore 
the nonlinearity is not confined to solely drive the slowly 
varying amplitude A. 



B. Transverse power flow reversing 

The extreme nonlinear regime can support propaga- 
tion of beams characterized by exotic properties such as 
the reversing of the electromagnetic power flow along 
the beam transverse profile [23. In order to discuss 
this phenomenon, consider a TM field describing a beam 
mainly propagating along the z-axis, i.e. a field of 



the form H y (x,z) 



A(x, z) with the requirement 



\d z A\ <C K\A\ (K being a wave- vector depending on 
the actual electromagnetic configuration). In this case 
Maxwell equations of Eqs. (fT0|) yield 



Hi, Dr 

v K * 



KE X + i 



dE z 
dx 



(16) 



where it has been assumed that d z H y ~ iKH y and 



d z E x ~ iKE x . Exploiting Eqs. ([T6|) . the ^-component of 
the time averaged Poynting vector S z = (l/2)Re(H y E*) 
can be expressed through the equivalent relations 



S z 

s z 



2K 



Re(D x E* x ) 



1 



-Re 



KE, r . 



dE z 
dx 



e: 



(17) 



If, for example, e > 0, u > 0, x < and 7 = 0, for an 
electromagnetic beam whose peak electric field strength 



is greater than \J\eJx\ (extreme nonlinear regime), it 
is evident from the first of Eqs.(fT2|) that D x and E x 
are antiparallel around the propagation axis (i.e. where 
|x(E • E*)| > e) and parallel elsewhere so that, from the 
first of Eqs. (fl7l) . S z is negative near the beam axis and 
positive elsewhere. In other words the beam is charac- 
terized by a power flow whose direction reverses its sign 
along the transverse profile and this is due to sign flip- 
ping of D x along the wave transverse profile while E x 
does not change its sign. Note that the sign flipping of 
D x corresponds to a sign flipping of H y (see the first of 
Eqs. p^]) ) so that, considering the second of Eqs.(fT6"|). this 
can happen without sign flipping of E x only if dE z /dx is 
not negligible with respect KE X . Therefore power flow 
reversing can take place only if the field has a transverse 
size comparable with 1/K. It is worth stressing that the 
discussed power flow reversing is very different from the 
effect that, in left handed metamaterials, the Poynting 
vector is antiparallel to the carrier wave vector which is 
a consequence of the fact that, in such media, e < 
and \i < (with n < 0). On the other hand, in our case, 
/i > and the sign of the power flow is not uniform being 
controlled by the field intensity through the nonlinearity. 



C. Linear and nonlinear parameters design 

In order to discuss the impact of the wide ranges of 
the parameters of Eqs. (|13p . made possible by linear and 
nonlinear design, on electromagnetic phenomenology, we 
consider the situation where the effective dielectric per- 
mittivity e and magnetic permeability p, are real, so 
that, after substituting the expression for the TM field 
of Eqs. pTj) into Eqs. fTU)) , eliminating the magnetic field 
and using Eqs. (fT2|) we get 

d 2 U z d 2 U x 



+ <r M [(1 + j)\u x \ 2 u x + (u x u: + ~fu:u z ) u z 



d 2 U v d 2 U, 



didC, di 2 

= a e a^U z + £r M [(1 + j)\U z \ 2 U z + (U Z U* + 7 U* Z U X ) U x ] , 

(18) 

where dimensionless variables and fields have been intro- 
duced according to 





(19) 



and fco — lu/c whereas a e = sign(ex) and c M = sign(//x) 
are the signs of the products ex and /^x> respectively. As 
opposed to the linear regime where electric field behav- 
ior solely depends on the sign of eu pj, from Eqs.(fl8]l 
we note that, in the present approach, the nonlinear dy- 
namics separately depends on the signs of ex and i- e - 
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the presence of the nonlinearity breaks the symmetry be- 
tween the roles played by the signs of e and /j, [H, [HJ • 
From Eqs.(fl"9f we note that e, /i and x scale the actual 
physical size and amplitude of the field so that, since the 
effective parameters can in principle be independently 
chosen (see Eqs. ([l3|) ). for each solution of Eqs. ([T8|) . a 
suitable layered medium can be designed in such a way 
that the actual electromagnetic field has prescribed geo- 
metrical size and intensity (see Eqs. dT9)) ). As an example, 
in the case of beam propagation, such an electromagnetic 
scaling freedom can allow to observe nonparaxial feature 
of a beam whose transverse width is much greater than 
the vacuum wavelength (if \efi\ << 1) or, on the con- 
trary, to observe the standard paraxial phenomenology 
for beams whose transverse width is much smaller than 
the vacuum wavelength (if |e/x| >> 1) and, remarkably, 
this can be done by avoiding any unfeasible requirement 
on the intensity. 

The parameter 7 plays a role fundamentally differ- 
ent since it cannot be generally removed from Eqs. (fT8)) 
with a field transformation. Setting U x = A a 
and U z = A z e % ^" (where A x , A z , (f> x and <p z 
real), the nonlinear terms of Eqs. (fT5)) , namely N x 
(1 + 7 )|t4| 2 C4 + {U x U; + -yU*U l ,)U a and N z = (1 
l)\U z \ 2 U z + (U Z U* + jU*U x ) U x , can be written as 



are 



N x 
N z 



(l + 7 )^+(l + 7 e- 2i ^-«)^ 



A x e^ 



(l + 7 e 21 ^"^)) A 2 X + (1 + 7) A 2 ] A z e 1 ^ , 



(20) 



from which we note that 7 is responsible for a nonlin- 
ear coupling between the two field components which is 
highly sensitive to the phase difference <fi x — <fi z . The 
case 7 = — 1 is particularly interesting since in Eqs. (l20[) 
the terms proportional to (1 + 7) vanish so that each 
field component is nonlinearly driven only by the other 
component and such coupling is fully sensitive the phase 
difference. Also intriguing is the regime I7I >> 1 since, 
two major opposite situations exist. If the phases of U x 
and U z are equal, from Eqs. (l2"0)l it is evident that the 
overall nonlinear polarization is proportional to (I + 7) 
so that, if I7I >> 1, the effective nonlinear response can 
be enhanced. If, on the other hand, the phase differ- 
ence between U x and U z is 7r/2, the overall nonlinear 
response of Eqs. ([2"0)l solely contains the terms (A 2 , + A 2 Z ) 
and 7 (A 2 — Al), so that, if I7I >> 1 (for realistic bounded 
fields) the compensation A 2 ~ A 2 has to occur in order 
to prevent the divergence of the term proportional to 7 
in Eqs.CEE]). 



IV. NONLINEAR GUIDED WAVES 

As explained in the above section, the nonlinear re- 
sponse of the proposed Kerr metamaterial is easy to man- 
age and, since the accesible ranges of its parameters are 
very wide, one can devise situations where the medium 



supports a genuinely novel nonlinear electrodynamical 
phenomenology. In order to probe the novel regime, we 
consider here a class of fields which are sufficiently simple 
to admit full analitical treatment and, at the same time, 
rigged with enough structure to show many of the novel 
effects. More specifically we consider nonlinear guided 
waves propagating along the z— axis of the form 



- JK 



u x (i), 



U z ((,0 = e^iu z (0, 



(21) 



where j3 is a real constant and the amplitudes u x and u z 
are real. Substituting the field of Eqs. (f2"Tj) into Eqs. (fT8|) 
we obtain 

-P—TT + P U x = We + U x , 

„du x d 2 u z . . . . 

P-jf -^ = ^[^ + *J ( 22 ) 

where 

*x(0 = (1 + 7)^(0 + (1-7)^(0, 

*»(0 = (l-7)«x(0 + (l+7)««(0- (23) 

Using Eqs. (f2"5|) , it is worth noting that, from the first of 
Eqs. (fT2l) . the vector D can be expressed as 



D x = e x 




{cr, + ^ x )u x e lK , 



D; = <t,\ \l - (a e + -f z )iu z e lK 

from which we note that 

JNL) _ 



(24) 



|e|sign(x) (cr e + * B ) 
4^) = k|sign( X ) (<r £ + *,) 



(25) 



act, for the fields of Eqs. (l21|) . as effective nonlinear di- 
electric permittivity. The quantities of Eqs.(|23p play a 
fundamental role in our discussion since it is evident that, 
if the conditions |* x | << 1 and \$f z \ « 1 do not hold 
along the profile of a nonlinear guided wave, the linear 
and nonlinear contribution in Eq.(|24p are comparable, 
i.e. the considered nonlinear guided wave belongs to the 
extreme nonlinear regime we have discussed in Sec lIII Al 
We consider solutions of Eqs. (|22|) with definite parity 
where u x and u z are spatially even (tt x (£) = Ux(~ 0) 
and odd (u z (£) — —u z (—£)), respectively and, as a con- 
sequence, we adopt the boundary conditions 



u x (Q) = u x0 , 
u x (+oo) = u xoc 



«,(0) 
u z (+oo) 



0, 

u 7A 



(26) 



It is worth stressing that, due to the feasible possibil- 
ity of arbitrary choosing the effective linear and non- 
linear parameters characterizing the effective nonlinear 
medium, we consider here solutions of Eqs. (l2"2l and (|2"6")l 
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TABLE I: Guided waves existence ranges of m^oo depending 
on 7, <r e and cr M . 





cr E = — 1 a^ = \ 


<X £ = — 1 cr M = —1 


7 > 




T 2-|_4 T -|_1 <- "io= 2 


-1< 7 < 






7 < -1 




7^ < «S < \ 



for <r e = ±1, (Tu = ±1 and, remarkably, for any real 7. In 



Appendix A we show that the system of Eqs. ([22|) is inte- 
grable and we derive the existence conditions character- 
izing the nonlinear guided waves satisfying Eqs. (|26p . i.e, 
for each possible combinations of cr e and tr M , we derive a 7 
dependent range of v 2 ^ (i.e. 1^(7) < < "max (7)) 
corresponding to nonlinear waves existence. The result- 
ing phenomenology is reported in Table I. For each ob- 
tained nonlinear guided wave the propagation constant 
/? and the asymptotical value u zoo are given by (see Ap- 
pendix A) 



P 2 = 



(<7 e + 2u 2 xoo ) 



1 + 7 

-cr e - (l-7)« 



■*x(+0o)], 



2 

x 00 



1 + 7 



(27) 



where the second expression for /3 is obtained from the 
first by exploiting the first of Eqs.(|23|). It is worth not- 
ing that, for each u xoo , these values of ft and u zoo are 
obtained by requiring that the nonlinear guided wave 
is asymptotically spatially uniform (i.e. by annulling 
the derivatives in Eqs. ([22|) for £ — > 00) and, exploiting 
Eqs.(|231). this implies that 



ANL), 



ei-"J(+oo) = ^(3\ 



0. 



(28) 



From a physical point of view, the second of Eqs. 
states that the considered nonlinear waves can propagate 
through the medium only if, asymptotically, the nonlin- 
ear term exactly balances the linear one to yield an over- 
all vanishing dielectric constant, i.e. 



* z (+oo) = 



(29) 



This proves that the nonlinear guided waves we are con- 
sidering always belong to the extreme nonlinear regime, 
i.e. they can not be observed in standard Kerr media. It 
is worth stressing that the sub-family of nonlinear guided 
waves with u zoo = do not require such an asymptotical 
compensation mechanisms since, for £ —¥ 00, the sec- 
ond of Eqs. dm vanishes together with u zoo . To sum up, 
an asymptotically spatially uniform wave can exist only 
if D z (oo) = (since d x H y = —iu>D z ) and this can be 

achieved by requiring either u zoo — or ei NL \+oo) = 0, 
the second situation being investigated in the present pa- 
per. 



Since the waves we are investigating always belong to 
the extreme nonlinear regime, power flow reversing dis- 
cussed in Sec lIIIBl is expected. For the field of Eqs. ([2~T|) . 
the time averaged Poynting vector S = (l/2)i?e[E x H*] 
is given by 



Sosign(/i) 



G e sign(x) , , , 2 - 
S = S — - — (a e + W x )u x e z 



(30) 



where So = \J (eo|e| 3 )/(4A*oHlxl 2 ) an d the second ex- 
pression is obtained from the first one by exploiting the 
first of Eqs. ([22|) to eliminate the derivative of u z and the 
the first of Eqs.([23]l. Note that Eqs.([30]l coincide with 
Eqs. pTl) (with the identification K = y/\efi\ko(3), the 
difference lying in the fact that, for the nonlinear guided 
waves we are considering, Eqs. ([5(H) are exact due to the 
waves propagation invariance. 



A. Nonlinear guided waves for cr E = — 1 and 07, = 1 



and 



= 1 is very intrigu- 



The situation a e = — 1 
ing since sign(e/i) = a e a^ = — 1 and therefore, in the 
absence of nonlinearity, no propagation can occur since 
the medium can support only evanescent waves. On the 
contrary, we have shown that, in this situation (see first 
column of Table I), the nonlinear Kerr metamaterial can 
support propagating nonlinear guided waves (i.e. with 
real propagation constant j3). In order to grasp this fact, 
we note, from the first of Eqs. ([2~T|) . for a t = — 1 and 
a,, = 1, we obtain 



/3 2 = -l + * a (- 



(31) 



from which it is evident that, in the absence of the nonlin- 
earity, /3 2 = — 1, i.e. only evanescent waves exist. How- 

and the second of 



ever, combining the first of Eqs.([23f 
Eqs. ([27]). we obtain 



* x (+oo) = 1 



27 
1 + 7 



1 - 2ui 



(32) 



from which it is evident that, for — 1 < 7 < and 
u 2 Xao < 1/2 (see the first column of Table I), *a;(+oo) > 1 
so that 1 > and the nonlinear wave can propagate 
through the medium. From a physical point of view, the 
same result can be understood by regarding the nonlin- 
ear guided wave as a background infinite plane wave with 
a distortion around £ = 0. The background nonlinear 
plane wave can propagate through the medium since its 
amplitude is such that \E , a; (+oo) > 1 so that the effective 
nonlinear dielectric permittivity in the first of Eqs. (|25|) 



is such that \itx NL ^ = [— 1 + ^ , 2; (+oo)] > 0. In other- 
words, in the extreme nonlinear regime, the nonlinear po- 
larizability overcomes the linear dielectric contribution in 
such way that the sign of the effective overall dielectric 
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FIG. 2: Nonlinear guided waves transverse profile of u x (£) 
(subplot (a)) and m z (£) (subplot (b)) at different values of 
u xao for a e = —1, er M = 1, 7 = —1/2. 



(a) (c) 




FIG. 3: Profiles of u x and u z of various nonlinear guided 
waves with <r e = — 1 and ct m = 1 and with 7 spanning the 
range —0.8 < 7 < —0.1. All the waves are characterized by 
the same asymptotical value u xoo = vO-4- For each 7 two 
waves exist, the first being reported in panel (a) and (b) and 
the second in panel (c) and (d) 



permittivity is opposed to that of the linear dielectric 
permittivity and the waves are consistently not evanes- 
cent. 

In Fig. (2) we plot various profiles of u x (^) and u z (£) 
for 7 = -1/2 and for u 2 xoo < 1/2. In Fig. (3) we plot the 
profiles of u x (£) and u z (£) for different nonlinear guided 
waves, in the range —0.8 < 7 < —0.1, each characterized 
by the same asymptotical value u xoo = yOA. As ex- 
plained in Appendix A, more than one nonlinear guided 
wave (for a given 7) can be generally found for each u xoo 
and, in the situation of Fig. (3), there are specifically two 



< 




FIG. 4: Profile of the normalized Poynting vector 
z— component S z /So evaluated for the nonlinear guided waves 
of Fig. (3). Here sign(e) = —1, sign(/i) = 1 and sign(x) = 1. 
Panel (a) contains S z /So evaluated for the waves reported in 
panel (a) and (b) of Fig. (3). Panel (b) contains S z /So eval- 
uated for the waves reported in panel (c) and (d) of Fig. (3). 
In the inset of panel (a), the width A (root-mean-square de- 
viation) of the various S z /So reported in panel (a) is plotted 
as a function of 7. 



waves, the first being reported in panel (a) and (b) and 
the second in panel (c) and (d) . The most striking feature 
emerging from Fig. (3) is that, for a given u xoo the closer 
7 to —1, the sharper the profile of u x , i.e. 7 produces an 
hyper-focusing effect for this waves when it approaches 
the value 7 = — 1. As expected (see the discussion of 
Sec lIII Cj) the situation 7 = — 1 is peculiar displaying 
a phenomenology absent for other values of 7. From a 
physical point of view the hyper-focusing effect results 
from the combination of two different mechanisms. For 
the first one, we note that, since e[ NL \+oo) = (see 
the second of Eas.(|28p). if 7 is very close to —1, from 
the second of Eqs.J^H) and the second of Eqs. (f2"3"|) , we 
conclude that u ZQO is much greater than one (in agree- 
ment with the second of Eqs. (f2"Tj) evaluated for 7 ~ —1). 
Therefore, |u z (£)l is much greater than one everywhere 
apart a small region around £ = where u z is very small 
(since u z (0) — 0) and du z /d£, is very large. The sec- 
ond mechanism supporting hyper-focusing is based on 
the fact that, for 7 ~ —1 the field a;— component (u x ) 
is nonlinearly driven solely by u z (see the discussion of 
Sec lIII Cj) . Therefore, in the region where u z is very small 
and du z /d(; is very large, from the first of Eas.(|22p. it 
is evident that u x (t;) displays a very pronounced peak, 
i.e. u x is tightly squeezed by u z for 7 ~ —1. In Fig. (4) 
we report the z-component of the normalized Poynting 
vector S z /Sq (see Eqs. fBTij) ) evaluated for the waves re- 
ported in Fig. (3). Note that the above discussed hyper- 
focusing effectively corresponds to a tight energy local- 
ization around £ = for 7 close to —1. This effect is 
particularly evident from the inset of panel (a) of Fig. (4) 
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where we plot the peak width A of S z /So as a function 
of 7. The power flows reported in panel (b) of Fig. (4) 
clearly displays the transverse power flow reversing dis- 
cussed in Sec IIII Bl since there is a region around £ = 
where S z /Sq < whereas S z /Sq > elsewhere. Note 
that the power flows reported in panel (a) of Fig. (4) does 
not exhibit transverse power flow reversing and this can 
be easily understood considering the structures of the two 
waves reported in Fig. (3). The z— component of the first 
wave is such that du z /dt; < (see panel (b) of Fig. (3)) 
so that, from the first of Eq. ([3U|) . the two bell shaped 
contributions u x and —u x du z /dt; are both positive yield- 
ing the S z /Sq > 0. On the other hand, the second wave 
is such that du z /d£, > (see panel (d) of Fig. (3)) and 
therefore, in the first of Eq. ([30|) . the two contributions 
have different signs and S z /Sq can flip its sign along the 
transverse profile. 



B. Nonlinear guided waves for cr e = — 1 and c M = — 1 

For (j e — — 1 and er M = — 1, waves are allowed to propa- 
gate in the linear regime and the nonlinear guided waves 
can propagate if 



2 



l-sa+oo) >o, 



(33) 



(see the first of Eqs.(|27p). i.e. the background nonlin- 
ear plane wave can not produce a nonlinear dielectric 
response overcoming the linear part. As opposed to 
the case discussed in Sec lIV A] in the present situation 
there are two distinct families of nonlinear guided waves 
corresponding to the ranges 7 > and 7 < — I (see the 
second column of Table I). In Fig. 5 we report different 
nonlinear guided waves profiles together with their power 
flows for different values of 7 > and for u xoo — y/OA. 
From panel (a) of Fig. 5 we note that, in this regime, the 
profile of u x has a peak and a hole (around £ = 0) for 
small and large values of 7, respectively. This behavior 
is easily understood since for 7 = we have /3 = (from 
the first of Eas.(|27|l) and u x (0) = 1 > u xoo (from the 
first of Eqs. (j22|) evaluated at £ = and for a e = —1) 
whereas, for 7 — > +00 we have (3 — > y/2(l — 2u|^) and 
Ux(0) — < u xoo (from the limit 7 — > +00 of the first 
of Eas.(|22p for a e = — 1). As a consequence there must 
be a value of 7 for which u x is uniform. Requiring that 
Ux(0 — u xoo, Eas.([2"2|) can be casted in the form 



(1 + 7) (-u 2 zoo +u 2 z )u z 



(34) 



so that the solution of the second of Eqs. (|3~4"f fulfilling 
the boundary conditions of Eqs. (|!2o]) is 



«.(0 



~,tanh 



1 + 7 



(35) 




FIG. 5: Profiles of u x (a), u z (b) and S z /So (c) of various 
nonlinear guided waves with er e = — 1 and = — 1 and 7 > 
0. All the waves are characterized by the same asymptotical 
value u xao = v/OA 
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FIG. 6: Profiles of u x (a), u z (b) and S z /So (c) of various 
nonlinear guided waves with cr e = — 1 and = — 1 and 7 < 
— 1. All the waves are characterized by the same asymptotical 
value u xao = \/0.4. 
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or, in other words, the longitudinal component u z is an 
exact electromagnetic dark soliton. However, the ob- 
tained u z has to satisfy the full Maxwell system so that, 
substituting the held in Eq. ([55|) into the first of Eqs. (j3"4"|) . 
it is straightforward to prove that this is possible only if 



1 



7± 



^ U XOO 



1± x l-Aui 



(36) 



For 

u xoo — \/0.4 we obtain 7+ — 2 which is the value of 7 
at which u x is uniform in panel (a) of Fig. 5. Note that the 
situation where u x is uniform and u z is a dark soliton is 
possible only in the extreme nonlinear regime where the 
nonlinearity ^ z can compensate the linear part in the 
second of Eqs. (|22|) . From panel (c) of Fig. 5 it is evident 
that every considered waves exhibit the transverse power 
flow reversing discussed in Sec lIIIBI 

In Fig. 6 we report different nonlinear guided waves pro- 
files together with their power flows for different values 
of 7 < and for u xoo = \/0.4. Let us consider, in this 
case, the behavior of the nonlinear guided waves for large 
values of |7|. From panel (a) and (b) of Fig. 6 it is ev- 
ident that a kind of saturation occurs, i.e. u x and u z 
approaches their asymptotic profiles for 7 — > —00. This 
is consistent with the fact that, taking the limit 7 — > —00 
of Eqs. ([27|) we obtain 



(3 2 = 2(1- 
u 2 = u 2 

zoo xoo 



2ut 



(37) 



i.e. ft 2 and u 2 zoo asymptotically approaches two finite 
asymptotic values and it is relevant that u 2 ZQO — > u xoo . 
In order to obtain the asymptotic profiles of u x and u z , 
we note that, taking the limit 7 — > —00 of Eqs. (f2"2"|) 
and consistently assuming that the profiles remains ev- 
erywere finite, consistency requires that u 2 (£,) = u1(0 
so that, asymptotically, u. x (£) = ti 2 (f) for £ < and 
u x(0 = —u z (^) for £ > (since u x and u z are spatially 
even and odd, respectively). Exploiting this property, 
Eqs. (|2"2")l yields, for 7 — > —00, the equation 



d u z 



4(- 



' u loo + u l) 



(38) 



so that, the asymptotic nonlinear guided wave profiles 
are 



u x{i) = Mxootanh (V2u xoo £*j 
(V2u xoo ^ 



Mz(£) = -Wxootanh 



(39) 



Note that the asymptotic profile of u x is singular at £ = 
(i.e. it continuous but not diffcrcntiable) as a conse- 
quence of the unrealistic assumption 7 — > — 00. Even 
though the profiles of Eqs. (l5^)) are not physical, they 



show that both u x and u z asymptotically approaches 
the square of an electromagnetic dark solitons and they 
provide a very good description of the above mentioned 
saturation for 7 < 0. The profiles of u x and u z for 




% 1 



FIG. 7: Profiles of u x and u z (solid lines) of nonlinear guided 
waves for u IOO = y/OA and 7 = —2, —10, —500. The asymp- 
totical (7 — > —00) profiles of Eqs. (139p are also reported (dot- 
ted line). 



u XO o = VOT and 7 = —2,-10,-500 together with the 
asymptotical profiles of Eqs. (|3"9")) are plotted in Fig. (7) 
from which it is evident that, for 7 < —10, the asymp- 
totical profiles of Eqs. (|3"9"]) accurately describe the actual 
field profiles (apart a small region around £ = for u x ). 
We conclude that, in the extreme nonlinear regime, the 
compensation between ty x and ^> z and the linear terms, 
if I7I >> 1, forces u 2 and u 2 almost to coincide. 



V. CONCLUSIONS 

In conclusion, we have investigated the effective re- 
sponse of a nonlinear Kerr metamaterial obtained by ho- 
mogenizing a one dimensional layered periodic structure. 
The effective response formally coincides with that of a 
standard nonlinear Kerr medium with the important dif- 
ference that its parameters (both linear and nonlinear) 
can be independently designed and that they can assume 
even values not achievable in standard material. As a 
consequence we can choose the linear dielectric permittiv- 
ity to be much smaller than one thus allowing the obser- 
vation of the regime where the nonlinear polarization can 
not be regarded as a small perturbation (extreme nonlin- 
ear regime) . As a leading general phenomenon character- 
izing the extreme nonlinear regime we have discussed the 
transverse power flow reversing effect, i.e. the fact that 
the power flow can change its sign along the transverse 
beam profile. Combining the extreme nonlinear regime 
and the fact the effective Kerr response can be tailored 
in an unconventional way, we have discussed a number 
of novel phenomena exploiting a class of fields (nonlinear 
guided waves) admitting full analytical description. Ex- 
amples of such novel effects are the fact that the nonlinear 
waves can propagate even if the medium linear proper- 
ties (dielectric permittivity and magnetic permeability) 
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would forbid propagation, hyper-focusing induced by the 
phase difference between the field components (in the 
case 7 = —1) and extreme compensation between the 
field components if they are 7r/2 out of phase in the lim- 
iting situation 7 >> 1. 



to cast Eqs. ([22|) into the standard form of a first order 
system of differential equations. Differentiating the first 
of Eqs.(|22|) and substituting the obtained expression of 
d 2 u z /d^ 2 into the second of Eqs. (j2"2"]l . we obtain, after 
some algebra 



Appendix A: Nonlinear guided waves existence 

In order to derive the existence conditions of nonlinear 
waves satisfying Eqs. (|22|) and Eqs. (f26]) . it is convenient 



y dU Z 

j dux_ 



(/3 2 - a e a,j,)u x - er M [(1 + 7)^ + (1 - j)u 2 z ] u x , 

(3 2 a e + (1 - 7) {(2^ - (3 2 ) + 2a, [(1 + 7)^ + (1 - 7)^] } u 2 x + fi 2 {\ + 7)^ 
a e + [3(l + 7 K + (l- 7 ) u 2] 



(Al) 



which is a system of ordinary differential equations equiv- 
alent to Maxwell equations provided the relation 

£r e +[3(l + 7)u| + (l-7)«g ^0 (A2) 

holds along the whole profile u x (^), u z (£) [35] . The sys- 
tem of Eqs. (|Al[) can be fully analytically investigated 
since it admits the first integral 

F(u x ,u z ) = [(/3 2 - cr e o>)u 2 - 0^07, u 2 ] 



~ { (/3 2 - (r e a M ) - a, [(1 + 7 )u 2 + (1 - 7 ) u 2 ] } 2 u \ 

(A3) 



or, in other words, the relation 



-F(u x (O,u z (O) = 
d£ 



(A4) 



holds for any solution u x (£), u z (£) of Eqs. (|Al[) . Consider- 
ing the boundary conditions of Eqs.(|26p. since u x (£) and 
u z (£) have to asymptotically approach two constant val- 
ues, their first and second derivatives vanish for £ — > +00 
so that we require the right hand sides of Eqs. (|A1|) to van- 
ish at u x — u XO o and u z = u zoo . As a consequence we 
obtain 



P 2 = ^(v e +2u xoo ), 
1 + 7 

,2 _ - C 1 - l) u . 



2 

x 00 



1+7 



(A5) 



which are the relations expressing the propagation con- 
stant P and the asymptotical amplitude u zoo as functions 



of the asymptotical amplitude u xoo . In addition, exploit- 
ing the fact that the F(u x , u z ) is constant along the wave 
profile, the relation 



F(u xO ,0) = F{u 



XOO 1 "zoo ) 1 



(A6) 



in which j3 and u zoo have been eliminated using Eqs. (|A5p . 
is a cubic equation for u 2 Q which can be solved to yield 
the peak amplitude u x q as a function of the asymptotical 
amplitude u xoo . 

For each u xoo , the existence of the corresponding 
guided waves has to be assured by a number of require- 
ments. In first place /3 and u zoo have to be real so that, 
from Eqs. (|A5p , we obtain 



1- 
(1 



J^kl >o, 



l) U lc 



1 + 7 



> 0, 



(A7) 



which are necessary inequalities for guided waves exis- 
tence. Analogously, it is necessary that u x0 , obtained by 
Eq. (IA6l) . is real. On the other hand, the above bound- 
ary conditions for £ — > +00 imply that, after substitut- 
ing the expression of (3 2 of the first of Eqs. (|A5P into 
Eq. (|A3l) . the function F(u x ,u z ) has a stationary point 
at Koo,«zoo)- Therefore, since the curve (u x (£),u z (£)) 
of the plane (u x ,u z ) has to reach the point (u XQC , u zoo ), 
it has to be required that F(u Xl u z ) has a saddle point at 
(«iooi"zoo) and this leads to the necessary inequality 



7 2 + 4 7 + 1 
1 + 7 



7CT e 



*°° 1+7 
X [2(1 ~ 7)<oo + ^(3 - j)u 2 XOD + 1] > 



(A8) 



We conclude that the existence of the nonlinear guided 
waves we are considering is assured by the condition of 
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Eq. (IA2l) . together with the reality of u x q (obtained by (IA8|) . The fulfilment of all these requirements leads to 
solving Eq. (|A6p ) and the inequalities of Eqs. (|A7[) and the existence conditions reported in Table I. 
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